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Abstract 

In this paper we reformulate free field theory models defined on the rectangular D + l dimensional 
lattices as D + 1 evolution models. This evolution is in part a simple linear evolution on free 
( "creation" and "annihilation" ) operators. Formal eigenvectors of this linear evolution can be directly 
constructed, and them play the role of the "physical" creation and annihilation operators. These 
operators being completed by a "physical" vacuum vector give the spectrum of the evolution operator, 
as well as the trace of the evolution operator give a correct expression for the partition function. As 
an example, Bazhanov - Baxter's free bosonic model is considered. 
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1 Introduction 

In this paper we investigate the simplest case of integrable statistical mechanics and field theory models: 
the free models. 

To fix some definitions and terminology, recall in few words, for example, scholar scalar free field 
theory in four dimensions. The hamiltonian 

H (k, <t>) = \ J dx (n(x) 2 + (d^x)) 2 + m 2 cb(x) 2 ) , (1.1) 

gives the evolution in physical time. In the quantum field theory a matrix element of the evolution 
operator can be presented in the form of the path integral of a Gaussian exponent. Hamiltonian ( [iTl| ) 
can be diagonalized by the Fourier transformation 

If - dk 

4>(x) = J (o+(-fe)+o(fe)) exp(ikx)^= (1.2) 



and 



= ( 2 7r)3/ 2 / e X p(ikx)i^dk, (1.3) 



with the dispersion relation 



= k 2 + mf . (1.4) 



The hamiltonian becomes 

H = J uj a+(k) a(k) dk , [a(k) , a+ (k')} = 5{k-k'). (1.5) 
The n - particle eigenvectors of the hamiltonian and so of the evolution operator are the Fock states 

n 

11 a + {kj) \0>, H |0>= 0. (1.6) 

3=1 

General free models on regular lattices have the same feature as this scalar model: all them can be 
diagonalized by proper Fourier transformations jlj]. The differences between the scalar model and a 
general free theory on the rectangular lattice are: 

• on the lattice the main object is not a hamiltonian, but an one-step evolution operator in the 
discrete time. 



• dispersion relations arising in the lattice models are little bit more complicated then (1.4), and they 
give nontrivial partition functions of the lattice models endowed by a phase transition phenomena. 

In the well known realm of the two dimensional models one could mention the free fermionic model 
(F.F.M.) as such case. The partition function for F.F.M. can be represented in a form of a path - type 
integral over the Grassmanian variables. The partition function is a determinant of large enough matrix 
which can be block - diagonalized with a help of the Fourier transformation, in this way the well known 
expression for the partition function for finite lattice can be obtained. [] 

An advantage of a 1 + 1 evolution formulation for the free bosonic and fermionic models is that 
one may avoid the path - type integration, considering the evolution of the creation and annihilation 
operators. 

It is very useful also to deal not with the usual transfer matrices as the evolution operators, but with 
the transfer matrices appearing in the diagonal sections of 2D lattice. Such matrices U do not form a 
commutative family, but they commute with the set of usual transfer matrices T, this is provided by the 
Yang - Baxter equation. 

lr The same is valid for the 3D F.F.M. iff 
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Fig. 1. The configuration U ■ T. 



U(^) ■ T{p,q) = T{p,q) ■ . (1.7) 

q q 

The graphic configuration is shown in Fig. 1. 

Such formulation of a free model as 1 + 1 evolution system with the cross - type evolution operators 
can be generalized easily to the (D + 1) - dimensional case with [D + 1) - simplex equation as the 
condition of the existence of a commuting T - type operators family. Moreover, for D + 1 evolution 
system it is also easy to construct eigenvectors for U and so for T. Note, the transfer matrices U in 3-D 
were investigated originally in Ref. In this paper we deal with 2 + 1 evolution models. 

The paper is organized as follows. In section 2 an operator formulation of free models in terms of 
Gaussian operators is given. In section 3 we define a 2 + 1 evolution model on the cagome lattice and 
describe a formal diagonalization of it. Examples are given in section 4. Most interesting example is 
connected with Bazhanov - Baxter's free bosonic model ||, and so section 4 is devoted to it mainly. 

2 Formulation 



Let A be an algebra with a left A- module F. Let A as the linear space has a basis e a . Let Q(A) be a 
group of linear automorphisms of A: for any G G Q {A) 



i = G o e Q = ^2 9a, p ep , 



(2.1) 



where G is properly defined on F so as 



G o e d = G ■ e ■ G' 1 



(2.2) 
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Call these G the Gaussian operators on A. To exhibit the definition of G via the matrix g a ^ we will use 
the universal notation G = G(g). Also we will omit the indices a denoting the basis of A via e so that 



G(g) ° e = g ■ e . 
Consider now the tensor product of A copies of F; 



F® A = F®F< 



F . 



(2.3) 



(2.4) 



A times 

This linear space is the left module of the direct sum of A copies of A. Let 

e k = l<g)l<g>...<g)e<g>l<g>.... 

I 



The Gaussian operators from G(%2 Ak), 

\ 

G(g) o 



( 9i,i .9i,2 

52,1 92,2 
\ 9A.1 9&,2 



k — th place 

91, A 

92, A 

5A,A / 



(2.5) 











/ 





(2.6) 



we suppose to be properly defined on F® A . 

Note that the matrices g represent the antihomomorphism of the group of the Gaussian operators 
G(g): 

G(g) ■ G(g') = G{g' ■ g). (2.7) 

Suppose for some reason we consider the Gaussian operators of special type: operators which act 
nontrivially only on some part of F® A . Namely, let 



then denote 

an operator such that 

while for e*k with k G h 

G h o 



h = {h, fc 2 , k M } , 

Gh = Gki,ka,...,kM 
Gh o ej — ej , j , 



efei 



5fci,fcj, 



\ 9k M ,k\ ■ ■ ■ 9k M ,k M 



Such operators appear, for example, when one considers the D-simplex equation: 



(2.8) 

(2.9) 
(2.10) 

(2.11) 
(2.12) 



where each ft( ra ' consists on D indices, the whole number of the indices in (2.12) is — 2 — an< ^ 
denoting 

h = {fa,..., fcu} ^ fc„ = (&)„ , (2.13) 
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then 



Note, h 



(2.14) 



All operators Gh are defined in End F® A , ( 2.12 ) is the usual D-simplex equation for the operators 
acting in the tensor product. But (2.12) implies also the corresponding equation for g^: 



(D) 



■ g'h> ■ 9h = gh ■ g' w 



(D) 



(2.15) 

and the main difference between ( 2.120 and (2.15) is that g^-s act on the direct sum of linear spaces: 
gh € End Ak- Such kind of D-simplex-type equations was investigated in Refs. Q. {efc} is the 

k 

formal basis, and G(gh)-S are defined on the vectors 

A 

cj = -e k . (2.16) 



k=l 



Thus there defined the co-action of G(g) on the amplitudes x k : 

x k h-> x' k = {g T ■ x) k 

due to 



G(g) o uj = ^ (g T ■ x)l ■ e k 



(2.17) 



(2.18) 



These simple facts allows one to formulate and solve eigenvalue and eigenvector problem of a special 
type. This will be done in the next section. 



3 2 + 1 evolution 

In this section we'll formulate some 2 + 1 evolution problem. The formulation is rather dimension 
independent, as well as the solution. We choose the dimension 2 + 1 because this is not so convenient as 
1 + 1 evolution problem and is not so cumbersome as an arbitrary D + 1 case. Also the examples we'll 
give belong to the 2 + 1 case. 

Suppose we are intending to investigate a model (statistical mechanics or field theory) defined on the 
cubic lattice: the cites of the lattice are 

ft = n x e x + n y e y + n z e 2 , (3.1) 

where e x ,e y ,e z are orthogonal normal basis in E3, n x ,n y ,n z € Zn, N is the size of the lattice, the 
periodical boundary conditions assumed. Instead of dealing with the usual layer - to - layer transfer 
matrices, i. e. section of the lattice between any n z = k and n z = k + 1, consider the section of the 31? 
lattice between n x + n y + n z = m and n z + n y + n z = m + 1. Call this section n x + n y + n z — m as Z m 
(do not mix this with the ring of residues) and define the evolution as the transition ... Z m 1— > Z m+ \ ... . 
Such evolution problem was suggested in Refs. ^j. 

Consider now the section, e. g. Z , in details. The two - dimensional oriented cagome - type lattice 
Z is shown in Fig. 2. Z consists of the oriented triples of the vertices, 

Z = { Triples [ V ait0j , V aulk , V 0j<lk ] } , (3.2) 

where the numbers are chosen so as i+j + k = 0, and V a ,(i stands for a vertex at the intersection of the 
lines a and (3 etc. 
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Oli+l Pj+l 




page. Fig. 2. A fragment of the lattice Z . 

Fix some notations for the sake of simplicity and shortness. First, introduce the 3D translations as 
the re-enumeration operators: 

t z ° [v ai ,i3 j ,v aulk ,Vf3 ] ^ k \ = \y ai+1 ,b j ,v ai+1 ^ k ,Vp jnk ] , 

T y o \y ai ,p j ,v olink ,Vp^ lk ] = [v ait b j+1 , v ai ,~f k , Vp j+Ulk \ , 

li"^,^,^] = [V ai , bj ,V ai ^ k+1 , Vp jilk+1 ] . (3.3) 
Pure 2D translations on the kagome lattice are 

Pa = T y T^ 1 , pfj = T z T^ 1 , p 7 = T z T' 1 , (3.4) 

so as 

P0 = Pa ■ P 1 ■ (3-5) 
Thus the position of any triple i, j, k : i + j + k = can be noted by 

P = Pa* ■ p\ (3.6) 
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because of 

P ° [V ao ,p , V aono , Vg , 70 ] = [V ait p } , V aulk , V Pj . lk } . (3.7) 

Note, by this notation we chose ao,/3oj7o to be special point in the frame of reference. 
Let the sizes of the lattice be N in each direction: 

p» = p% = p» = 1 , (3.8) 

so that the total number of the vertices is A = 3 N 2 . 

Assign now to each vertex V the algebra Ay with the basis eV as well as its module Fy ■ Then for a 
triple [V at p, Vq i7 , Vp n ] introduce a (Gaussian) operator R: 

R : F Va ® F Van ® F Vp „ h-» Fy 0i<J ® Fy a , 7 ® F Vp „ , (3.9) 

where it is supposed that in the right hand side the orientation of the vertices [V aj p, V a ^, Vp a \ is changed 
as it is shown in Fig. 3. 




page. Fig. 3. The action of R. 

To distinguish the configurations call the left hand side one as L.H.S. triple and the right hand side 
- as R.H.S. triple. Consider now the evolution of Z generated by applying R-s to all L.H.S. triples of 
Z . As the result we obtain Z\ formed by other L.H.S. triples: 

Z x = { L.H.S. Triples [ V ai3j , V Q „ 7fe , Vp jl7k : i + j + k = 1] } . (3.10) 

Applying R-s to new L.H.S. triples repeatedly one obtains 

Z -> Z x -> Z 2 -» Z 3 -» ... (3.11) 

where 

Z m = { L.H.S. Triples [ V a<fij , V aulk , Vp jnk : i + j + k = m ] } . (3.12) 
For U m : Z m Z m+ \ one may write the following expression : 

Um = n r v „ v v a ■ ( 3 ' 13 ) 

-L.L v ai,f3ji v ai,jki v fijHk 

i+j-\-k—ra 
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Introduce shorter notations for ey and Fy assigned to the vertices V ail ^, V aulk and Vp jt7k . Namely, 
mark everything for V a< p by x, for V an - by y, and Vp n - by z, and the position of a triple [V at p, V an , Vp n ] 
- by its p as in Eq. (|3.6|): 

V au0j = pi ■ p}) o V ao ,p = V^pijfp) <-> 

e,^), ^(p^), (3.14) 

^,7* = Pa^ ' P% ° V ao , yo = Vyfe^) «-» 

eyip-%), F y (p-%), (3.15) 

? z ( P p k p-3), F z (p- k p-r), (3.16) 
With these notations the action of i?, Eq. ( |3.9D , looks like 

R : F x (p) ® F y (p) (g> F z (p) ^ F x (t xP ) ® F v (t vP ) ® F z (t zP ) . (3.17) 
This situation is shown in Fig. 4, analogous to Fig. 3. 




page. Fig. 4. The action of R. 

Let now all operators R coincide: 

r v*( p ),v v (p'),v4p") = R(ex(p),e y (p'),e z (p")) (3.18) 

with the same operator function R. Then all U m differ only by some re-enumeration of the vertices. 
Choose the re-enumeration operator r = t x : 

T O V X ( P ) = V x (p) , 

T O Vy{p) = Vyip-'p) , 

r o V z (p) = V.(pp l p) . (3.19) 

Then, obviously, 

U m — T m o Uq — T m ■ Uq • T~ m . (3.20) 
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The m - layers evolution operator [/( m ) : Z n Z m , 

(7 (m) = U m -i ■ U m - 2 ■ ... • Ui ■ U = r m ■ (r- 1 • U ) m . (3.21) 
Thus the universal one step evolution operator 

U = t- 1 ■ U : Z a ^ Z (3.22) 

arises. 

Now we are going to diagonalize U. 

Let further R-s (and so U) are indeed Gaussian operators. Investigate the action of R and U on the 
linear space ^ Ak of 

e x (p) e y (p) , e z (p) ■ (3.23) 

Let the nontrivial part of 

r eQ(A x @A y @A z ) , (3.24) 
R = R(r), is given by 3 x 3 block - matrix 





( e x 








)■( 




v 1 





' x.x i 1 x,y •) 1 x.z 



R(r) o e y = r v>x , r ViV , r y>z • e y . (3.25) 



e. 



The space Ax (p) © A/ (p) © A (p) consists of the vectors 

w = E % ■ + ' + ' , (3-26) 

p 

where the amplitudes x p ,y p , z p may be chosen so that 

U o lo = Xlu . (3.27) 



Obviously, 
Hence 



and 



= k% = 1. Thus 



p Q • (7 - U ■ Pa , p p ■ U = U ■ p/3 . (3.28) 

p a o uj = ^ %p ■ e x (p a p) + ... 
p 

= }2 ■ + - = ka w ' ( 3 - 29 ) 

V 

pp o co — kp uj , (3.30) 

^V^p^p = & a [x,y,z] p , = kp[x,y,z\ p . (3.31) 



Consider now the co-action of U on the amplitudes (see Fig. 5): with the eigenvalue A 



' Z,X 1 



—*± \ -*l I —*± — #/ 

f J A = ^p • + jf-i p • r ^ + f J-V • • (3.32) 
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U 




Up z -i U -i 

Pa P "Pa, P 



page. Fig. 5. Transformation of the amplitudes. 

Take into account the translation symmetry, then the finite eigenvector and eigenvalue problem arises: 

- \ \ T / 



r y,x, r y,y > 

r z,yj Tz,z *^/3 J \ 

Solving this problem, one obtains the eigenvalues of U as the roots of the characteristic polynom (the 
dispersion relation) 

X^Afc^Afc^lr) = & \i = \i(k a ,kp) , i = l,...,3dinU, (3.34) 

where 

// Ar \ \ 
x(Ai,A 2 ,A 3 |r) = det A 2 - r . (3.35) 
VV A 3 J J 

Let further 

\n>e F® A (3.36) 

be the vacuum vector for U, i. e. 

U ■ |fi > = \Q> (3.37) 

and \Q > is the cyclic vector in F® A . Suppose e are free elements of A, then the set of eigenvectors of 
U is the following set of the Fock - type vectors 

|m,...,n A >= w£ ■ ... ■ wj£ • \n> (3.38) 

with the eigenvalues 

U ■ K,...,n A > = A? 1 ... A A a |m,...,n A > • (3.39) 

Then formally 

Trace U M = J] , (3.40) 



Afc-\ r y . z • y p = . (3.33) 
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so for the partition function per cite zo(r) the bulk free energy ko(r) = logzo( r ) hi the thermodynamic 
limit is given by 



1 N ~ 1 

M»') = - L im iv.M^oo jpjj log Il ^ ~ X f ( fc "'M) 



J /• /• /-27T 



(27T) 



k a ,k b =0 j 

d(f>i dfo dfo log X (e #1 , , |r) . (3.41) 



In general vectors ( 3.3 are not all independent, in this case one should take into account the diffe rence 
between "the number of the degrees of freedom" and the dimension of "phase space" multiplying ( 3.41 ) 
by an appropriate constant fc l— * k — rj k . More precisely, for a given representation of on F one 
should extract from the set of the eigenvectors lo\ — u>i(k a ,kp) the operators, cyclic with respect to 
|fi >, and sum the power series of their eigenvalues. 



Thus for a Gaussian local evolution operators one can easily construct the eigenvectors as well as the free 
energy for the thermodynamic limit. Note that the local integrability conditions (D - simplex equations) 
we have not used at all, because the diagonalization of our U and the existence of a set of commutative 
operators are not linked at all. Nevertheless, the cases of the integrable models are more interesting 
and we'll deal with the evolution operators U constructed with the help of R-s solving the Tetrahedron 
equation. 



4 Functional R — operator 

The interesting example of the Gaussian operators connected with integrable models can be obtained 
from the functional R - operators solving the functional tetrahedron equation (F.T.E.) Q. 
Let a functional operator R: 

#1,2,3 [x 1 ,X 2 ,X 3 ] = [x[,x' 2 ,x' 3 ] , (4.1) 

where 

x\ = ri(x 1 ,x 2 ,x 3 ) , x' 2 = r 2 (xi,x 2 ,x 3 ) , x' 3 = r 3 {x\,x 2 ,x 3 ) , (4.2) 

solves F.T.E. Examples of such R-s one can find in Ref. Mention now a very useful aspect of the 
functional maps. Namely, let Xk E Sk- Then actually R can be upgraded to R giving an isomorphism 
on a bundle of Si £g> S 2 ® S 3 . Partially, on the tangent bundle, 

Ro [xi, x 2 , x 3 \dxi, dx 2 , dx 3 ] = x 2 , x' 3 \dx[, dx' 2l dx' 3 ] , (4-3) 

where 

dx 'i -E^^> ( 4 - 4 ) 

so as we can introduce the Gaussian operator for the differentials on S: 

r = R ■ R- 1 , (4.5) 

where, using the symbol a, instead of dxi, 

dx' 

r o en = a- = ^ r i,k{x) afc , n t k(x) = R^ 1 o ^-i . (4.6) 
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Thus for each functional operator R : S® 3 i— > S® 3 there appears the gaussian action G(r) : dS i— > 
d$ for a still commutative algebra dS. can be raised to nontrivial dS © 

r o aj = £ 4 (r- 1 )*,* - (4-7) 

with the conserving Weyl algebra [ai, a^] = 

Because of i? solves F.T.E., G(r) solves a local T.E. in the direct sum of dS © dS © .... Thus the 
previous considerations give the solution of the eigenvalue and eigenfunction problem for the 2 + 1 
evolution governed by the properly quantized R. 

Among other models one could mention several free bosonic or fermionic models, namely Bazhanov 
- Stroganov's free fermionic model and Bazhanov - Baxter's free bosonic model. 



4.1 Bazhanov-Baxter's free bosonic model 

Recall in few words the form of R - operator connected with the complex of the descendants of the 
Zamolodchikov - Bazhanov - Baxter model @, ||, ||, |], [ToJ, |ll| . Let 

u ■ v — q v ■ u , w = —q^ 1 / 2 u ■ v . (4-8) 

Then for the quantum dilogarithm function J^, |l(| 

oo 

i){ X ) = JJ (1 - q 1/2+n X) (4.9) 
n=0 

over u,v,w there exists the Pentagon equation 

ip(v) ■ tp(u) = ip{u) ■ ip(w) ■ ip(v) , (4-10) 
and as the sequence of it the Tetrahedron equation for 

Rx,y,z = V'O" 1 W z )ll}{w~ 1 U z )n x ^^ll){w x U~ 1 Y 1 IpiWyW' 1 )- 1 , (4.11) 

@, where 

Rx,y,z ° {v x ,V y ,V z } = {v x ,V x V z ,V~ 1 V y } , 
Rx,y, Z O {u Xl U yi U z } = {u x WyW^ X 1 V~ 1 U z ,V x U y } , 

Hx,y,Z {W X ,Wy,W Z } = {W X WyW~ 1 ,W Zl Wy} . (4-12) 

Profs. Bazhanov and Baxter in 1992 established the correspondence between the pair u, v and the 
Gaussian exponents of the canonical Weil pair a, a + || []. First, show as the Gaussian action like (4.5) 



This correspondence on the language of w — functions on the Fermat curves is given by 

— » i da , w ab — + cxp(-iab) , 



_A(v) 

3=1 

i 



i \ TT &(y) i v 2 
■J(v,a)= w(v,a) = exp{--- a ) 

-1- J- 1 — VU)3 2 1 — V 



*(o) = (-) V / 2 - *(a) = exp(--a 2 ) , 

A(v) N = l-v N -» A(v) = 1 - v . (4.13) 

With this correspondence all the key relations in the Zamolodchikov - Bazhanov — Baxter model conserve, W - type model 
is well defined and its partition function is proportional to the partition function of usual finite - states ZBB model. 
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arises in the language of the quantum dilogarithms. First, easily, 

tp(u) ° v — v ■ (1 — q 1 / 2 u)~ x 

rp(v) o u = u ■ (1 - q~ 1/2 v) (4.14) 
In the limit when q 1 / 2 —> 1 i/j(u) and tp( v ) become functional operators ipf(u) and ipf(v) and 

= i« ' iA u )f , i>( v )f = L v ■ ip(v)f , (4.15) 

where for a + = — and a = — , 

It D 

a • a+ = a+ ■ a + A . (4.16) 
'+' is not h.c, q = exp(A), the Gaussian operators 



<• f „+ ) « ( i ) - ( ^ ) • ( : + ) 



where from the automorphism condition [a',a' + ] = A it follows that aS — f3j = 1, arise. Eqs. ( |Q]| ) 
with 

L„ = ex P(^ a 2 ) and i« = ex P(^ a+2 ) ( 4 - 18 ) 

give 

-<£° 2 > • ( ; ) - ( i : ) • ( :+ ) • <« 9 » 
-<£« +2 > • ( ; + ) - ( j i ) ■ ( :+ ) ■ <«■» 

These operators should be completed by the diagonal one, 

«" )^( <1+ " r ; (1+r >M )• ,421) 

where the sign :: of the normal ordering is used: formally 

: exp(^a+ • a) : = V a+" • a™ . (4.22) 

A i n! 

These three operators give the Gaussian triangular expansion of CP(2). Mention one more operator, 
r , I a » / 1 — r — r \ f a 



ex P (- (a+a+ no^ ; + j = ^ ; ;+ j. (4.23) 

The exponential forms (these Gaussian exponents define our terminology) are rather formal, because 
their explicit forms depend actually on a representation, and the following integral form is more useful 
sometimes then the series decomposition: 

exp(|^a 2 ) = J d(exp(^C 2 + ip<Za) (4.24) 

etc. 

There are two remarkable relations for these Gaussian operators: 
first, Kashaev's relation [ p"2| : 

exp(|a 2 ) • exp(|a+ 2 ) • exp(^a 2 ) = exp(|a+ 2 ) • exp(|a 2 ) • exp(|-a+ 2 ) , (4.25) 
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where the electric network transformation arises 



xy 



x + z — X 2 xy z 1 
y = x + z — X 2 xy z , 

, v z 

z = . 

x + z — X^xyz 

Second, Bazhanov's "restricted star - triangle" for the bosonic 3D model ||: 

exp(^a 2 ) • exp(^a+ 2 ) = exp(^-a+ 2 ) • exp(^-( a + + a) 2 ) ■ exp(^-a 2 ) 



"2X 



2X 



where 



*2A 



1-P 

q = q 



2A' 



<2X 



P = P 



1-pq' 
1-9 



1-pq 



r = pq 



(4.26) 



(4.27) 



(4.28) 



This is the usual pentagon relation connected with eq. ( 4.10 ). Although we can't extract expa directly 
from the quantum dilogarithm, the identification of the spectral parameters of exp-s with the functional 
transformation generated by the quantum dilogarithm ip(v) when g 1 / 2 — > 1 and v i— ► v exp(a), following 
from eq. (4.15), gives the correspondence 



Func. Part tpf(u)] ■ exp(— 



1 v 
2Xl-v 



(4.29) 



Give just an example a representation of eq. (4.27). Let A be pure imaginary so that the pair a and 
a + is the pair coordinate - momentum. Take the representation when 



so that 



and 



Then 



d 

a ■ \x > = (A x — — ) [a; > 



a + ■ \x > 



0_ 

dx 



\x > 



< x\ ■ a = < x\ , < x\ ■ a + = < x\ (Xx — -^-) 
ox ox 



<y\x> = exp(- x 2 ) 5(x - y) 



1 = J \x > exp(——x 2 )dx < x\ 



< x\exp(-^a 2 )\y > 
< x\exp(^-a +2 )\y > 



cxp(-— (x - yf 
2p 



exp(-— (a:-?/) + - 



2 V) 

X 



r X 
< z|exp(— ( a + +a) 2 )\y > = exp(- (1 + r) y 2 ) S(x - y) 



(4.30) 

(4.31) 
(4.32) 



(4.33) 



In this representation there appears exactly Bazhanov - Baxter's bosonic restricted star - triangle 
relation. 
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Return now to R - matrix (4.11). We'll identify 



v — > v exp(a — a + ) , u — > it exp(a + ) , w — > to exp(a) 
We'll deal with the endomorphisms of the direct sum of several pairs of 



at 



(4.34) 



(4.35) 



The equivalence between ip and expa 2 is given by (4.29). 

As it was described in the previous sections, for the R - matrix being treated as the Gaussian operator 
R = G(r), we are interested in the matrix r and the evolution governed by this r. Thus for n Xj2/jZ 



1 , 7T„ 







and 



Restore all other exp-s: 



1 
1 




-1 1 



Rx 



y.z 



exp(;rT-K - azf) ■ exp( — (a x - at) 2 ) 



'2X 

P3 



T x ,y,z ■ ex P(-^( a 2; - a z 



2A 

CX P(-^( a y - a z) ) 



where from eq. ( 4.11 ) it follows that 

Pi(l - Pi) = Pi{l - p 3 )(l - Pi) ■ 

Two matrix functions arise: 



and 



Then 



P 



> cr(p) 



a(p) 



and 



</(p) = 



and 



/ 1 














o \ 




p 


1 











-p 










1 






















1 










p 











1 


-p 




V o 














1 J 




/ i 













\ 





1 



















1 
















P 


1 




-p 
















1 


V o 





-p 







p 1 


/ 




P3) 




,y,z 




cr(p 2 ) • ( 



(4.36) 
(4.37) 



(4.38) 

(4.39) 
(4.40) 
(4.41) 

(4.42) 



(4.43) 



(4.44) 
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Thus the eigenvector problem for the bosonic model reduces to the six - dimensional eigenvector problem. 



Let amplitudes for eq. (4.35) be 

xl = {x p , x*) etc. (4.45) 
(do not mix p as a point on the lattice with p as a parameter) . Then the amplitudes of the eigenvectors 



are given by eq. ( 3.33| ), where the final expression for r (4.44) is 



/ 1-Pl 





Q2 





-Q2 


Pi 


\ 





i-Qi 


P3 


Qi 


-P3 


o 




Pi 





I-Q2 





O2 


-Pi 




-Qa 


Pi 





1-P2 





Q 3 







-Qi 


P 3 


Qi 


1-P3 







V Q 3 


-Pi 





P 2 





1-Q3 


) 


Pi = P2 , 


P2 = 


1 - P3P4 , 


P3 = 


i-p 3 (i 


-pi) , 




Ql = P3 , 


Q2 = 


1 - PlP2 , 


Q3 = 


i-p 2 (i 


-Pa) ■ 





(4.46) 



(4.47) 



The characteristic polynom (3.35) for r - matrix (4.44) does not factorize except some special limits, 
and the direct calculations for the integral (3.41) give the expression 



where 



and Roger's dilogarithm 



Zq(pi,P2,P3,Pa) = Bfa) B(p 2 ) B(p 3 ) B(p A ) , 

B(p) = exp -e(p) (Li(p) - Li(l/2)) , 
e(p) = sign Im(p) , 

p ^ log x log 1 — x 



Li(p) 



1 

2 ,„ 



1 - x 



(4.48) 

(4.49) 
(4.50) 

(4.51) 



These calculations are cumbersome, but they can be simplified significantly in the limit when pi =Pa = 0. 
In this case the characteristic polynom factorizes and there appeared the integral definition of B: 



B(p) 



1 



exp - 



d<t>i #2 d<t>3 ■ 



(2^) 3 j j jo 

log(l - p(e* + e lW2+ ^ 3) ) - (1 - p){e^ 2 + e 1 ^ 1 ^) + e ^+^+^) 



(4.52) 



This integral can be taken easily and gives ( 4.49| ). Match the expression for zq with Baxter's results for 
Zamolodchikov model z 2 jlq] and Bazhanov's result for the bosonic model zb py: 



(4.53) 



The origin of the difference between zq and Bazhanov's zb is that up to the proper representation our 
system of the eigenvectors is overdefined twice (see the remark at the end of Secti on 3) . 



Few words concerning the eigenvectors of U. First, the characteristic polynom ( 3.35 ) for r ( 4.46 ) has 
the following property: 

X(x,y,z\r) = x 2 y 2 z 2 x{x~ l ,z~ l \r) . (4.54) 



Thus for the isotropic eigenvectors x = y = z six roots of x = divide into two subsets: 

■ 1 ^ - J ^ -A3 1 } 



{Ai, A 2 , A 3 } and {A 4 = A x \ A 5 = A 2 1 , A 6 



(4.55) 
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shuch that | Ax 2,3] < 1 (note, |Afe| = 1 only if all p-s are real). This property remains for arbitrary 
quasimomenta k a ,kp: 

|Al,2,3(fea, Ml < 1 and ^4,,5fi( k a,kp) = Al, 2 ,3 (K 1 , k^ 1 )^ 1 . (4.56) 

In order to obtain the convergency of Trace U , suppose the eigenvectors for the linear evolution 
operator LO\.2,z{k a , kp) to be the creation operators, so the vacuum vector is to be defined 

^4,5,6(^/3) |fi> = . (4.57) 

Due to Eq. ( 4.56| ) the structure of the bulk free energy formula (3.41) for ko does not change, but the 
^ 

extra multiplier ^ appears: k — — ko. Recall, in the case then Pi,P2>P3>P4 are real, all A-s belong to 

the unit circle, hence when one changes the signs of Imp, the structure of the physical vacuum changes, 
so B(p) depends on the sign of Imp so drastically. 

Note, the derivation of the partition function via the three dimensional integral was suggested origi- 
nally for the free bosonic model by V. V. Bazhanov fll4|. 



5 Summary 

In this paper we have formulated the evolution problem for D + 1 dimensional free models. The evo- 
lution U - type operators act linearly on the algebra of the creation and annihilation operators, and 
so due to the locality of U the eigenvector problem can be solved directly up to a finite eigenvector 
problem. Appeared physical creation and annihilation operators depending on D momenta diagonalizc 
the evolution operator. In the case of integrable models, when there exists a set of T - type transfer 
matrices commuting with U, the eigenvectors diagonalize these T as well. The existence of the set T for 
the given U is provided by the D + 1 simplex equation Fl 
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